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Using ab initio density functional theory and quantum transport calculations based on nonequi- 
librium Green's function formalism we study structural, electronic, and transport properties of 
hydrogen-terminated short graphene nanoribbons (graphene flakes) and their functionalization with 
vanadium atoms. Rectangular graphene flakes are stable, having geometric and electronic structures 
quite similar to that of extended graphene nanoribbons. We show that a spin-polarized current can 
be produced by pure, hydrogenated rectangular graphene flakes by exploiting the spatially-separated 
edge states of the flake using asymmetric, non-magnetic contacts. Functionalization of the graphene 
flake with magnetic adatoms such as vanadium also leads to spin-polarized currents even with sym- 
metric contacts. We observe and discuss sharp discontinuities in the transmission spectra which 
arise from Fano resonances of localized states in the flake. 

PACS numbers: 73.63.-b, 72.25.-b, 75.75.+a 



I. INTRODUCTION 

In the last two decades various forms of carbon nanos- 
tructures, namely buckyballs, carbon nanotubes, and 
lately graphene have attracted a great deal of interest due 
to their novel fundamental properties and possible appli- 
cations in electronics. With the developments in prepa- 
ration and synthesis techniques carbon-based nanostruc- 
tures have emerged as one of the most promising mate- 
rials for non-silicon electronics. 

The fabrication of graphene sheets-L and observation of 
their unusual properties such as a half-integer quantum- 
Hall effect have attracted much interest in electronic 
transport properties of this type of two dimensional 
graphitic materials. Observed gapless energy spectrum 
and high mobility electron transportSii^ii^ are the most 
remarkable features of graphene. It was shown by tight- 
binding calculations considering the vr bands that in the 
electronic energy dispersion of graphene, energy is lin- 
early dependent on the wave vector around the Fermi 
level^ which makes it a unique material. 

In recent experimental studies graphene nanoribbons 
(GNRs) with narrow widths (10-70 nm) have been 
realizedi^ii Li et al report producing ultra narrow rib- 
bons with widths down to a few nanometers.— In addition 
to high carrier mobilities that are higher than those in 
commercial silicon wafers, existence of width-dependent 
energy band gaps makes the graphene nanoribbons a po- 
tentially useful structure for various applications. The 
width dependence of the band gap and transport prop- 
erties in quasi-one dimensional narrow GNRs have been 
studied theoretically^>iaai>i2a3>i4a5a6^ 

Graphene nanoribbons having r zigzag rows (r- 
ZGNRs) are predicted as semiconductors having a nar- 
rowing band gap with the increasing width of the rib- 
bon. Armchair-edged ribbons (AGNRs) are also semi- 
conducting with direct band gaps.^'' Density functional 
theory (DFT) calculations predict a high density of 
states around the valence and conduction band edges of 



ZGNRs, which derive from the states that are localized at 
the edges of the ribbon, and lead to non-zero magnetic 
moments on the carbon atoms. The existence of such 
peculiar edge-states in graphitic materials had been dis- 
cussed in some earlier workSi ^^'^^i^°i^^i^^'^'^ It has been 
established that the ground state of a ZGNR is anti- 
ferromagnetic with ferromagnetically ordered spin polar- 
ization along the edgesj2 Due to the fact that the states 
near the Fermi level are derived from the edge states and 
their linear combinations, external fields have more sig- 
nificant effects on the edge states. The nature of the 
interactions between the magnetized edges of ZGNR was 
studied, concluding that the main contribution to the 
local moments comes from dangling bonds, and their in- 
teraction is determined by tails of the edge-localized tt 
statesi^ Due to the edge effects, graphene ribbons show 
different one dimensional transport properties from those 
of carbon nanotubes. 

Since GNRs have long spin-correlation lengths and 
good ballistic transport characteristics they can be con- 
sidered as a promising active material of spintronic 
devices i^i^i^^i^ In particular ZGNRs, known to have 
large spin polarizations at the opposite edges of the rib- 
bon, may be utilized to create spin-dependent effects such 
as spin polarized currents without the need of ferromag- 
netic electrodes or other magnetic entities. 

Hydrogen termination of the edge atoms by forming 
strong a bonds would be important for the structural and 
electronic stability of the graphene ribbon. Both first- 
principles and tight binding calculations showed that the 
termination of edges with hydrogen atoms removes the 
electronic states related to the dangling bonds How- 
ever, there are no qualitative changes in the electronic 
structure and the magnetic order of the ZGNRs with hy- 
drogen atom terminationji^ except for a narrowing of the 
band gap as we calculated. 

Modification of electronic structure by impurities, 
adatoms and external fields is another potential of 
graphene and graphitic structures for applications in 
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nanoelectronics. Effects on electronic properties and 
magnetic beliavior of graphene by the adsorption of 
foreign atoms has been considered in some previous 
^Qj.]jg^29j30|3i|32|33|34j35j36 Sizeable gap opening by hydro- 
gen adsorption to Stone- Wales defect sites of graphene 
has been reported.— In the case of substitutional boron 
atom^^''^'^ many of the electronic properties have been 
studied and it is suggested that GNRs may be used as 
spin filter devices4^ Based on the spin polarized ground 
state of ZGNRs their possible application as a spin- valve 
device was proposed by using tight binding j^^i^^ and k • p 
calculations.^ - It is also predicted by using first-principles 
calculations that the electric current flowing on the rib- 
bon can be made completely spin polarized under in- 
planc homogeneous electric fields.— 

In this study we consider electronic transport prop- 
erties of short graphene nanoribbons (graphene flakes), 
and their functionalization with vanadium atoms. Elec- 
tronic and magnetic properties of various types of 
nanometer-sized graphene flakes have been reported 
before^^i^ii^i^^i^iSi^'^ We restrict our considera- 
tions, however, to rectangular flakes only, which are finite 
segments of perfect GNRs. In particular we calculate 
spin-dependent transmission spectra of the flakes when 
the electrodes make partial contacts along the zigzag 
rows of the flakes. The geometrical asymmetry of the 
contacts lead to a polarization in the spin states of trans- 
mitted electrons. Depending on the contact geometry 
and the electrode thickness the spin polarization of the 
transmission shows various forms around the Fermi level, 
including perfect (100%) polarization at certain energy 
ranges. Then we perform first-principles calculations to 
investigate the optimized geometry and electronic prop- 
erties of the flakes with adsorbate vanadium atoms. We 
discuss the binding sites for single and double adsorbate 
atoms on the graphene flake. By using the relaxed geom- 
etry of adsorbate-flake system, we calculate their conduc- 
tance spectrum by using model metallic electrodes. The 
adsorbate vanadium atoms introduce additional states 
around the Fermi level modifying the conductance spec- 
trum of the flake by breaking the spin symmetry of the 
conduction electrodes even when the electrodes make uni- 
form contacts with the flake. We also find traces of Fano 
resonances in the conductance spectrum of the graphene 
flakes arising from the coupling of the extended states 
of the quantum channel with the localized states of the 
flake, which acts like a scattering quantum dot. A local 
density of states (LDOS) analysis show that the discon- 
tinuity features with sharp peaks and dips in the con- 
ductance spectra are associated with perfectly localized 
states of the structure in the electrode-flake-electrode ge- 
ometry. 



II. CALCULATION METHODS 

We performed first-principles total energy calcu- 
lations to obtain electronic structure and equilib- 
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FIG. 1: (Color online) Geometry and spin-dependent charge 
density of the graphene fiake cut from 4-ZGNR. The edges 
are saturated with hydrogen atoms. Green (dark) and yel- 
low (light) regions denote the local majority spin-type of the 
charge density. Possible adsorption sites of adatoms are la- 
belled as El, H1-H4. Electrodes (carbon chain or gold bar) 
are attached to the flake along the directions A, Z1-Z4. 



rium geometries of rectangular graphene flakes, with 
and without adatoms, based on the pseudopoten- 
tial density-functional theoryi^ The spin-dependent 
exchange-correlation potential is approximated within 
the generalized gradient approximation^^ (GGA). The 
software package Atomistix ToolKit^ (ATK), which em- 
ploys local numerical basis orbitals and noncquilibrium 
Green's function formalism to calculate quantum conduc- 
tance in electrode-device-electrode geometry, has been 
used in all total energy and transport calculations. 

Geometry optimizations of graphene flakes have been 
done by relaxing all atomic positions in supercell geome- 
tries with sufhcient vacuum spaces (10 A minimum) 
to prevent the interactions with periodic images of the 
structure. The criteria of convergence for total energy 
and Hellman-Feynman forces were 10""* eV and 0.005 
eV/A, respectively. The electrostatic potentials were de- 
termined on a real-space grid with a mesh cutoff energy of 
150 Ry. Doublc-zeta double-polarized basis sets of local 
numerical orbitals were employed which gives consistent 
results with previous calculations^* that used plane- wave 
basis sets. In cases of extended nanoribbon calculations 
the Brillouin zone has been sampled with (51,1,1) points 
within the Monkhorst-Pack k-point sampling scheme. ^- 

The geometric structure of the graphene flake consid- 
ered in this study is shown in Fig. [TJ It is a finite seg- 
ment of 4-ZGNR structure having 4 zigzag rows of carbon 
atoms along its axis. The sides along the cut direction 
have armchair shape. All edge carbon atoms have been 
saturated with hydrogen atoms for a better structural 
and electronic stability of the flake. Possible adsorption 
sites of adatoms have been denoted in the figure. The 
magnetic moments of the atoms and the structure as a 
whole have been calculated by using Mulliken analysis. 

The spin-dependent transport properties of the 
graphene flake were calculated based on non-equilibrium 
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Green's function formalism as implemented in ATK. Bal- 
listic conductance spectrum of the flake is given within 
Landauer-Biittiker formalism as 

GiE)=GoJ2T.{E) (1) 

a 

where Gq = e'^/h is the conductance quantum and T„ 
{a = up, down) are the spin-dependent transmission 
spectra. For each spin state a the transmission prob- 
abilities are calculated as 

T,{E) - Tr[rLG''rRG'']„, (2) 

where C and are the retarded and advanced Green's 
functions, and Tn are the contact broadening func- 
tions associated with the left and right electrodes, respec- 
tively. The broadening functions are anti-Hermitian com- 
ponents of the self-energy terms of the contacts, S l and 
Tin, which depend on the electrode's surface Greens func- 
tion and the contact-molecule (graphene flake) bonding: 
^L,R = T'{'Sl,r — jj)/2. The spatially-separated spin- 
states of the graphene flake are anticipated to break the 
spin symmetry in these conductance calculations when 
contacts to the electrodes are not symmetric at the arm- 
chair edges. 

In order to exploit spin-dependent characteristics of 
the graphene flake the structure needs to be partially 
contacted with thin metallic electrodes. We used lin- 
ear carbon atomic chains as model electrodes. Carbon 
atomic chains which are known to be metallio^ are ex- 
pected to make reasonably good contact with the flake. 
The possible contact sites of the electrodes are denoted 
in Fig. [TJ Semi-infinite carbon atomic chains have been 
rigidly attached to the flake at one or more sites along the 
A, Z1-Z4 directions after removing the hydrogen atoms 
at the contact sites. Optimum electrode-flake distances 
were obtained by relaxation of the supercell consisting 
of the flake and 4-buffer carbon atoms per each carbon 
chain attached on both sides of the flakes. 

As an alternative to carbon-chain electrodes we also 
considered gold electrodes to test the robustness of our 
conductance calculations. The gold electrodes were de- 
vised by taking a finite cross-section bar of atoms oriented 
in the [001] direction of fee crystal. The semi-infinite gold 
bars have a unit cell of nine atoms in two layers of four 
and five atoms. We have added single gold atoms to 
the ends of the electrodes to get atomically sharp tips. 
Contact regions of the electrodes to the flake were geo- 
metrically optimized. 

III. RESULTS AND CONCLUSIONS 

For the generation of spin polarized currents using a 
graphene flake, formation of spin-ordered edge-localized 
states found along the zigzag edges is the key mechanism. 
Therefore, we first calculate the ground state magnetic 
ordering in ZGNR structures and compare it to that of 



the flake. In their ground state, all the carbon atoms of 
ZGNRs have antiparallel magnetic moments with their 
nearest neighbors, and their magnitudes rapidly decrease 
as we go away from the edges. For instance, accord- 
ing to the MuUiken population analysis moments of the 
atoms in the edge zigzag row of a 4-ZGNR alternate 
between I.21/Zb and -O.lO/ie per atom. In the second 
zigzag row the corresponding values are O.IQ/^b and - 
0.07/iB. Relatively high magnetic moments are found 
on the atoms that are under-coordinated. Upon satura- 
tion of these carbons with hydrogen atoms the moments 
on the edge carbon atoms are reduced but the qualita- 
tive properties of the magnetic ordering of the ZGNRs 
persisted. For comparison, in the hydrogen-saturated 4- 
ZGNR the edge-row carbon atoms have alternating mo- 
ments of 0.26^B and -0.06/Xb. In the second row the se- 
quence of moments reduce to 0.07(1^ and -0.04/iB- The 
hydrogen atoms too attain minute moments of about 
0.01/iB, aligned antiparallel to those of the carbon atoms 
they are attached to. The spin-polarized edge structure 
of hydrogenated ZGNRs persists for wider ribbons, even 
with a slight increase in the moments of edge atoms. Hy- 
drogenization also decreases the energy band gap of the 
r-ZGNRs. For instance, we calculated that the band gap 
values for r =3, 4, 5 reduce from 1.25, 1.15, 1.00, to 0.70, 
0.65, 0.60 in units of eV, respectively. 

Geometric and electronic structures of finite segments 
of ZGNRs, namely rectangular graphene flakes, are quite 
similar to those of extended ribbons. Fig. [1] shows the 
particular graphene flake (C52II20) considered for the cal- 
culations in this study. Additional armchair edges of the 
flakes are also saturated with hydrogen atoms. Spin- 
unrestricted ground state of the flake corresponds to the 
antiferromagnetic spin coupling between the zigzag edges 
of the flake. The total energy of this state is 0.31 eV 
lower than that of the ground state of a spin-restricted 
calculation, which is consistent with the calculations on 
extended ZGNRs.^'i We find that away from the arm- 
chair edges of the flake the magnetic moments of the 
zigzag-edge carbon atoms are close to that of infinite 
ZGNR structure, but tend to decrease toward the arm- 
chair edges. For instance, the moments of the six car- 
bon atoms located at the upmost zigzag edge of the flake 
shown in Fig. [T] are 0.09, 0.21, 0.29, 0.29, 0.21, 0.09 in 
units of fi^. The magnetic moment profile of the carbon 
atoms rapidly converge to the ribbon values for longer 
flakes. 

We calculate the spin-dependent ballistic conductance 
of the graphene flake when contacted with metallic elec- 
trodes. Since zero bias conductance and transmission 
probability spectra are related by a simple multiplicative 
factor in Eq. 1, we frequently use these terms as equiva- 
lent quantities in the following discussions. By exploiting 
the spatial separation of spin states at the opposite zigzag 
edges of the flake we show that spin-polarized currents 
can be achieved without using ferromagnetic electrodes 
and magnetic entities. By making partial contacts to 
the armchair sides of the flake the spin symmetry in the 
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FIG. 2: (Color online) Spin-dependent transmission spectrum 
of the graphene flake without (top panels) and with (mid- 
dle panels) hydrogen saturation of the edge atoms. The car- 
bon linear chain electrodes are connected to the Zl sites (see 
Fig. [1} . The spectrum when the carbon chain electrodes are 
replaced by gold leads are shown in the bottom panels. Fermi 
level is set to zero. 



coupling of the flake's electronic states to the electrodes 
is broken. Consequently the transmission spectrum of 
the electrons through the flake becomes dependent on 
their spin state. For a clear demonstration of this effect 
stripped off complications that may arise from electrode- 
graphene interactions we chose to use simple, model elec- 
trodes. Carbon linear atomic chains are appealing in that 
respect, being metallic, making naturally good contact 
with graphene, and having atomically sharp tips. 

In Fig. [5] we plot equilibrium transmission spectra of 
the bare and hydrogenated graphene flakes with single 
carbon atomic chain electrodes asymmetrically contacted 
to one of the edge zigzag rows of the flakes. For com- 
parison, the spectrum of hydrogenated flake with more 
realistic gold electrodes is also presented. Presence of 
dangling bonds in the bare flake is reflected in the trans- 
mission spectrum as several sharp peaks about 0.7-1.0 
eV above the Fermi level. These peaks are removed af- 
ter hydrogenization of the flake, and are replaced by a 
single broad peak degenerate for both spin channels. We 
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FIG. 3: (Color online) Spin-dependent transmission spectra 
of the hydrogen-saturated graphene flake when the electrodes 
make partial contact at the Zl and Z2 sites (top panels), and 
when the contacts are uniform (bottom panels). Fermi level 
is set to zero. 



should note that within the chosen convention of labelling 
spin states the electrodes are contacted at the flake edge 
where up-spin states were localized. Therefore, the up- 
spin states around the Fermi level couple more effectively 
with the electrodes leading to formation of a broader 
peak for up-spin channel in the differential conductance 
spectrum, and consequently the spin polarization of the 
current for small bias voltages is expected to reflect the 
flake's local majority-spin type of the contact side. With 
the gold electrodes, the basic features in the transmission 
spectrum remain unchanged except for the transmission 
peaks being narrower. The narrower peaks are due to 
weaker coupling of carbon atoms to gold electrodes. 

The effect of electrode thickness and its contact lo- 
cation on the transmission spectrum is further demon- 
strated in Fig. [31 When the carbon electrodes make a 
wider contact with the hydrogenated flake at Zl and Z2 
sites the spin asymmetry in the conductance spectrum 
is blurred with broader and overlapping peaks, but still 
the up-spin channel has a larger contribution to the con- 
ductance of the flake at energies around the Fermi level. 
In case of uniformly contacting electrodes (Zl to Z4) the 
transmission spectrum is identical for both spin channels 
as expected. Similarly, conductance of the flake when 
the electrodes are contacted along the armchair direction 
(contact site A in Fig. [1]) is spin symmetric (not shown). 

Another mechanism to break spin symmetry in the 
conductance of graphene flakes is to introduce magnetic 
impurities or adatoms. Here we present a systematic 
study of adsorption of vanadium atoms on graphene flake. 
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TABLE I: Calculated properties of vanadium adatoms on the 
graphene flake. See Fig.[T]for the notation of adsorption sites. 
The magnetic order in the flake before the adsorption is such 
that the magnetic moments of the carbon atoms at the top 
edge are up {fj.c > 0). In cases of double adsorbate atoms their 
locations relative to the flake have been denoted by top (t) 
or bottom (b). Magnetic moments of the vanadium adatoms 
(/iv) and the total moments of the structures ^tat are in units 
of Bohr magneton (/ib). £b is the binding energy per vana- 
dium adatom, and do-v is the equilibrium distance of the 
adatoms to the plane of graphene flake. 



Adsorption site 


Mv 


Mtot 


£b (eV) 


dc-v (A) 


El (/iv < 0) 


-3.8 


-4.9 


1.59 


1.94 


El (/iv > 0) 


3.8 


3.0 


1.58 


1.94 


H2 (/iv < 0) 


-3.9 


-4.3 


1.01 


1.86 


H2 (/(V > 0) 


3.9 


4.4 


1.02 


1.87 


H1-H2 (t-t) 


-2.5, 2.0 


0.0 


2.09 


1.79, 1.71 


H1-H3 (t-t) 


3.6, -3.6 


0.0 


1.42 


1.77, 1.77 


H1-H4 (t-t) 


-1.9, 1.8 


0.0 


2.41 


1.74, 1.74 


H2-H2 (t-b) 


3.5, -3.5 


0.0 


1.12 


1.75, 1.75 


El-El (t-b) 


3.7, -3.6 


0.9 


1.28 


1.98, 1.96 



Single atoms of several elements are adsorbed on carbon 
nanotube surfaces with relatively high binding energies 
and magnetic moments<^ Similar properties can be ex- 
pected for their adsorption on graphene. Vanadium is a 
suitable element to consider due to its high binding en- 
ergy and spin polarization. The ground state electronic 
configuration of vanadium is [Ar]4s^3d^, and it is a chem- 
ically active element with a magnetic moment of 3/iB. 

We investigated the adsorption of single or double 
vanadium atoms on the graphene flake. Possible adsorp- 
tion sites of vanadium atoms on the graphene flake are 
shown in Fig. [1] When there are two vanadium adatoms 
their location relative to the flake are denoted by top (t) 
or bottom (b). Our results are summarized in Table I. 
Binding energies per adatom are calculated by using the 
expression = {E[G]+nE[V]-E[G+nY])/n. Here, n is 
the number of adsorbed V atoms, -^[G] and -^[V] are the 
total energies of the graphene flake (C52H20) and a free 
vanadium atom, respectively. Total energy of the relaxed 
structure when foreign atoms are adsorbed on graphene 
is denoted by E[G+nV]. Like most of the 3d transition 
metal elements (Sc, Ti, Cr, Mn, Fe, Co) vanadium atoms 
are adsorbed on hollow sites of two dimensional graphene 
sheet. Due to the edge effects, however, on the graphene 
flake vanadium atoms are not adsorbed on the HI, H3, 
or H4 sites but rather they migrate to the El site after 
relaxation. We found that vanadium atom can be ad- 
sorbed also on the H2 site, but the binding energy is less 
than the one at the El site. H2-site adsorption is also off- 
center; equilibrium position of the vanadium atom tends 
toward one of the hexagon sides. 

Adsorbed vanadium atoms modify the spin-relaxed 




FIG. 4: (Color online) Majority spin-type of local charge den- 
sity of the flake with single vanadium adatoms. The upper 
zigzag edge of the flake was spin-up polarized before the ad- 
sorption. The vanadium atom has a down magnetic moment 
in (a) and (c) and an up moment in (b) and (d). The lowest 
energy conflguration for a single vanadium atom is shown in 
(a). 

charge density of the flake, however the planer geometric 
structure of the graphene flake is quite robust to vana- 
dium adsorption, the displacements of even the nearest 
neighbor carbon atoms are less than 0.05 after geometry 
relaxation. As shown in Table I depending on the initial 
spin configurations of the flake and vanadium we iden- 
tify two distinct magnetic states with close energies for 
each adsorption site of El and H2. In its ground state 
a vanadium atom adsorbed at the El site has a mag- 
netic moment antiparallel with the moments of the edge 
carbon atoms as depicted in Fig. Ua). 

Binding of two adsorbate vanadium atoms on the same 
side of the flake is always stronger than that of being on 
different sides. For both single-sided and double-sided 
cases, two vanadium adatom configurations with antipar- 
allel magnetic moment have larger binding energies than 
the cases with parallel magnetic moments. 

By means of magnetic dipol-dipol interactions stronger 
bindings occur between the edge carbon atoms and vana- 
dium adatom. Considerable change in the spin charge 
density of flake with vanadium adsorbate atoms [see 
Fig. |4] also modifies the conductance spectrum of the 
structure. In Fig. [5] it is seen that the adsorption of a 
vanadium atom at the El site with spin-down magnetic 
moment enhances the down-spin transmission channel 
around the Fermi level. Spin polarization of the con- 
ductance spectrum is present even when the electrodes 
make uniform contacts with the flake. 

In Fig. [S] we present an alternative view of the spin 
asymmetry in the conductance spectra by the plots of 
energy-dependent relative spin-polarization of the trans- 
mission probabilities, (Tup - T^^q^^) / {Tup + T^own) ^ cal- 
culated for four different cases. Partial contacts of the 
electrodes and presence of magnetic adatoms both pro- 
duce spin polarizations in the spectra. When the flake 
is contacted from the Zl sites either with carbon chain 
or gold electrodes the polarization is almost 100% for an 
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FIG. 5: (Color online) Spin-dependent transmission spectra 
of the hydrogen-saturated graphene flake with an adsorbed 
vanadium atom at the El site in its ground state. Top panels: 
The electrodes make partial contacts at the Zl sites. Bottom 
panels: When the contacts are uniform. 



energy window of ~0.8 eV around the Fermi level. We 
also see that while uniformly contacted flake does not 
exhibit spin polarization, it can be generated by a single 
vanadium adatom. 

Before we conclude, we would like to discuss several 
discontinuity features found in the transmission spectra, 
a clear example of which can be seen in top panels of 
Fig. [3) These are in the form of a pulsative variation with 
very narrow dips and adjacent peaks in the transmission 
function, reminiscent of a resonant effect. Similar reso- 
nances have been recently reported in the transmission 
spectra of carbon nanotubes with a single iron or vana- 
dium adsorbed atom on themi^ The underlying mech- 
anism of these variations in the transmission function is 
the Fano resonances of localized states of the structure 
with extended band states of the electrodes with match- 
ing symmetry. The origin of the localized states con- 
tributing to the resonances can be either the localized d 
states of the adatoms as in Ref. 61 and in the vanadium 
adatom parts of the present study, or more generally they 
can be the localized states of a quantum dot coupled to 
the continuous states of a quantum channel,-? For in- 
stance, in formation of the resonances found in Fig.[3]the 
graphene flake acts like a scattering quantum dot coupled 
to the quantum channel of carbon-chain electrodes. 

We infer about the character of these resonances by 
calculating local density of states (LDOS) of the struc- 
ture at resonance and off-resonance energies of the trans- 
mission spectra. In Fig. [7] we present LDOS isosurface 
plots for up-spin transmission spectra of three particular 
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FIG. 6: (Color online) Spin polarization in the transmission 
spectra of the graphene flake. The type and contact geome- 
tries of the electrodes are given as insets. In the lowest panel 
there is a single vanadium adatom in its ground state. 



contact geometries. LDOS are perfectly confined at the 
flake region only at resonance energies. At off-resonance 
energies we always find contribution of extended states 
of the metallic electrodes. 

In summary, we have performed an ab initio study 
of spin-resolved conductance spectra of rectangular 
graphene flakes. We have shown that spin-polarized cur- 
rents can be obtained by making the electrodes contact 
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the flake partially along the zigzag-edge directions of the 
flake, or by introducing transition-metal adatoms such 
as vanadium to alter the symmetry in spin-dependent 
scattering rates of the transmitted electrons. In the par- 
ticular model system considered in this study (a C52H20 
flake cut from 4-ZGNR) a 100% spin-polarization in the 
conductance spectrum was achieved for both model car- 
bon atomic chain and more realistic gold bar electrodes. 
Traces of Fano resonances have been found in the trans- 
mission functions. We believe that the qualitative fea- 
tures of producing spin-polarized currents in graphene 
flakes can be verifled experimentally with flakes of larger 
sizes and with realistic electrodes. An understanding and 
realization of spin-dependent effects in graphitic struc- 
tures are essential toward the development of graphene- 
based spintronic devices. 
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